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Pion structure functions are considered in the framework of QCD sum rules at the 
next to leading order. The 2-loop diargam, representing the gluon exchange between 
the struck quark and spectator, is calculated for space-like pion momentum p 2 < 0. 
lO ■ After Borel transformation we obtain a leading twist correction to the pion structure 

function, which does not violate Bjorken scaling. Partial results can be used in the 
■ calculation of complete set of a s corrections to the pion structure functions, which may 

. serve as an intitial condition in the evolution equation. Probabilistic interpretation of 

£^ [ the structure functions is discussed. 

in 
o 

g ; 1 Introduction 

The nucleon structure functions, namely, the parts of the 2-point electromagnetic current 
■ ! correlator 

jg : w ^ = hJ d ' zel9Z ^M m ( z )jt m mN( P )) (i) 

are measured in the process of deep inelastic scattering. In the Bjorken limit they represent 
the sum of the probabilities fi(x) to find a parton i carrying the part x = Q 2 /2(pq) of the 
nucleon momentum p. This long time known statement is based on the Wilson Operator 
Product Expansion, which leads to the so-called twist expansion (for a review see [T]): 

where Q 2 = —q 2 , dimensionless constants C originate from local operators Ml ... Mj of dimen- 
sion d and spin J, averaged by the nucleon state: 

(N(p)|O m ... MJ |JV(p)) =p„.. . VilJ M d - J - 2 C + ... (3) 

The most important are the lowest twist biquark operators O ~ qj^q- Introducing the light- 
cone wave functions, one comes to a probabilistic interpretation of the nucleon structure 
functions at high Q 2 . 
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The leading order contribution of the 2-quark operator is commongly visualized by the 
following graph: 

7* 'Wvvv fxAA/W 



N (4) 

The lower lines stand for the soft nucleon constituents (spectator), but the upper one shows 
hard (struck) quark, interacting with the virtual photon. QCD interaction has several topo- 
logically distinct diagrams: 



a) evolution 



b) initial state c) final state 

interaction with spectator (5) 



The interaction of the struck quark with gluon field (jSK) leads to renormalization of the 
operators (jHJ). It yields the InQ 2 dependence of the structure functions, governed by the 
evolution equation. Various anomalous dimensions have been calculated at the 2-loop order 
|21 , and the 3- loop results were published recently [3]. The diagrams (03), (JSt) with spectator 
involved display the contribution of the quark-quark-gluon operator. Consequently, they are 
either a gauge artefact jl] or higher twist corrections [E] . They do not admit any probabilistic 
intepretation. 

However, recently a statement appeared in literature [0], that initial or final state in- 
teraction with spectator may give rise to a leading twist corrections, so the propabilistic 
interpretation of the structure functions should be reconsidered. In [6J it was demonstrated 
by a model calculation of the difractive dissociation cross section. A source of these correc- 
tions is the P-exponent between the quark fields, necessary to keep the gauge invariance of 
the structure functions. Although its argument vanishes in the light cone gauge A + = 0, 
it may generate a nonzero contribution due to additional pole l/k + of the gluon propaga- 
tor here, as the authors of [6] demonstrated. For a more discussion of the P-exponent role 
see [7|. The initial and final state interactions were also shown jH] to give a leading-twist 
contribution in Drell-Yan process and semi-inclusive deep inelastic scattering. 

Effect of the spectator interaction can be explicitly calculated in the framework of the 
sum rule approach. Here the nucleon state is approximated by some current with appropriate 
quantum numbers, acting on the vacuum. Then the problem is reduced to the calculation of 
the 4-point current correlator in vacuum. It has poles at p 2 equal to the masses of physical 
states (nucleons) and discontinuity at the continuum. By means of dispersion relation it 
can be analitically continuated to space-like region p 2 < 0, where it can be calculated as 
perturbative QCD series. Borel transformation by p 2 suppresses the contibution of the higher 
states and continuum and one obtains explicit expression for the structure functions in the 
intermediate x region, where the sum rules are supposed to work. Such a program was 
performed at the leading order in jH] for spin 1/2 nucleons. 

The calculation of the a s corrections to the nucleon structure functions will involve 3-loop 
graphs. It is simpler, however, to begin from the pion case, since it has one loop less. The 
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probabilistic interpetation, anomalous dimensions and the statement of [H| do not rely on 
the hadron structure, so they are applicable to the pion as well. The subject of this paper 
is the calculation of the final state interaction in the pion. In the next section we briefly 
review the derivation of the pion structure functions at the leading order, performed in [TUj . 
and then present original results. Technical details are dropped to Appendix. 



2 Leading order pion structure functions 

One considers the 4-point current correlator: 

IW(p,g) = - J rfa;^^e^ + ™-^(0|Tjt(x)j- m ( 2 /)jr(0)j5/ 3 (^)|0) (6) 

where j$p = U'jpjsd is the pion axial current, = X]q e g<?7/i<? electromagnetic one. If 
momentum p is close to the pion mass shell p 2 ~ m 2 , one may insert the intermediate pion 
state, obtaining: 

= ( /_ P "J )2 / dyeto(K(p)\T 3 ™{y)jr(0)\*(p)) (?) 

where f n is the pion decay constant (0|j 5( g(0)|7r(p)) = if n p@. Authors of [TU] took the 
incoming and outgoing pions with different momenta p\ ^ p\ (but (piq) = (j>2<?))- However, 
it is not necessary for our purposes: at the leading order results are equal. Moreover, the 
2-loop calculation, performed in the next section, seems untractable in this case. Correlator 
((7j) has discontinuity by the variable s = (p + q) 2 : 

PaPB 

Disc s il^ a/3 = n MI/a/ a(s + iO) - ^^(s - iO) = ^j^TZ^2y W ^ v ^ 

The tensor consists of the pion structure functions. 

In the quark basis the correlator © is given by the box diagram at the leading order: 



7* '^^a/vv 
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(9) 



and the diagram with m, d quarks exchanged. It is easy to show, that all other diagrams 
(unlike the photon-photon scattering) are not essential in the Bjorken limit. 
Calculation of the diagram (JHl) gives the following imaginary part: 

3 e 2 x (i _ x ) / 4 X 2 a a \ 

Disc,n^(p, q) = - — p2 p a pfi I - e„ v + — P^P V J + h.t. (10) 

where = ^- , P^ = p M - ■ 
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Letters "h.t." stand for the higher twist contributions, i.e. ~ In (— p 2 ), 1, p 2 etc, small with 
respect to the leading term 1/p 2 in the Bjorken limit p 2 — > 0. 

In the QCD sum rule approach one applies Borel operator to the variable — p 2 in order 
to suppress the contribution of the higher states and continuum. Then the following sum 
rule can be obtained: 



B M 2Disc fl Il 



-ml/M 2 



PaPpW, 



M 4 

^z(l-ar) p a pp 



Ax 2 - - \ 
e^u + ~gi p ^ ) + h -t- ( n ) 



ttM 2 

The Borel mass should be taken JU] 

M 2 = Sirf 2 (12) 

Then one extracts the w-quark distribution in the pion: 

f u (x) = Qx(l-x) (13) 

In fact, the sum rules are not valid for x close to and 1, so exact value of M 2 as well as 
the normalization integral 



dxf u (x) 



(14) 



are modified by various effects. We will not pursue this point further, see ^U] for details. 
The sum rule (jllj) demonstrates the following properties of the parton model: 

• Callan-Gross relation F 2 = 2xF\. 

• Bjorken scaling F(x,Q 2 ) = F(x). 

• Probabilistic interpretation of the function (|T3~j) is supported by the normalization 
integral (JT3J). 

The Bjorken scaling does not survive the a s -corrections. In particular, the leading order 
evolution equation gives the InQ 2 term: 



A/, 



a: 



evol 



2 — lnQ*[l-x + 4(1- x)x\n(l- x) - 2{l-2x)x\o.x 
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(15) 



Vacuum condensates also have their contributions, see ^H] . In this paper we will pay atten- 
tion to perturbative calculation of the interaction with spectator. 



3 Gluon exchange with spectator 



In this section we will calculate the diagram 



p 

o 
o 
o 
o 



P 



(16) 



describing the interaction between the struck quark and spectator. To avoid extra disconti- 
nuities, the pion momentum is taken to be space-like p 2 = — P 2 < 0. 

We will consider its contribution to the tensor structure p^PvPaPfii i- e - to the function 
F 2 . It it separated out by multiplying the diagram (JTfiJl on the projector: 

o 4 r - - .... 

d(d _ 2)A 4 [3A 2 e^e a0) + Q{d+l)q 2 Ae { ^P a P p) + (d + l)(d + 3)q 4 P^P u P a Pfs 

where A = (pq) 2 — p 2 q 2 . One can put the space-time dimention d — 4 from the beginning 
since, as explained below, only finite integrals will be important. Then the calculation is 
reduced to evaluation of some scalar integrals. 

Such 2-loop scalar integrals were considered in the calculation of the a 2 anomalous di- 
mensions of the structure functions 0. There it was sufficient to put p 2 = from the very 
beginning, since all IR poles cancel in the operator renormalization. In the sum rule ap- 
proach one may not accept such simplification, since the leading order integral already 
behaves as 1/p 2 . In dimensional regularization all the divergences have the pole l/(d — 4), 
regardless of their structure (UV or IR). For instance, it seems impossible to distinguish the 
terms 1/p 2 and ln(— p 2 ) in such way. The first terms are interesting for us, while the later 
ones represent a higher twist correction. Consequently, we have to keep p 2 finite. 

The calculation is simpler in the Feynman gauge, since there are only 2 graphs of the 
leading twist: 



k + q I + q 





k — p I — p 



(17) 



Indeed, consider the left graph of eq (fTTj) . Let us choose a frame p = (0, P, 0j_) and take a 
part of the phase space with the integration momenta k,l ~ p. The volume of this region is 
J d 4 kdH ~ p s . Each propagator, depicted by the solid line on the graph, brings the factor 
1/p 2 here, while the dashed line brings 1/q 2 . Consequently the left graph has the order 
P 8 1 (p 2 ) 5 (<? 2 ) 2 — IAp 2 *? 2 )- The second graph of the leading twist, the right one in eq (JT7J), is 
obtained by removing single dash line. (If both dash lines are removed, the integral becomes 
real, since no intermediate state can be inserted.) But if any solid line is removed, the twist 
increases and such diagram can be ignored. 

However, the most complicated 2-box graph of (|T7|) turns out to be real: 

DlSC " l k ,l k 2 I 2 (k - I) 2 (k - p) 2 (I - p) 2 [k + q) 2 (I + q) 2 = ° (18) 

We do not know a simple way to prove this result, so we put (rather cumbersome) details 
of the calculation in Appendix. In the evaluation of (fTo^) one encounters 2 integrals of the 
box-triangle leading twist topology: 

Ji = Disc, / TTTrTr ^ 1 - TTTl — 7 : NO (19) 



■h 



Disc s 
Disc, 



k>l k 2 I 2 (k - I) 2 (k - p) 2 (I - p) 2 (I + q) 2 

2{kq) +q 2 

, { k 2 I 2 (k - I) 2 (k - p) 2 (I - p) 2 (I + q) 2 



(20) 
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Exact results for these integrals (arbitrary p 2 , q 2 , s) are available in Appendix (|A21|) . (|A22|) . 

The contribution of the diagram (fTtj)) to the structure p^PvPaPp of the correlator (jHJ) 
discontinuity is: 



2 11 7re 2 u a s ix{\ - x) [ - (1 - 2x) J x + x Q~ 2 J 2 ] + h.t. 



(21) 



Substituting it into the sum rule (fTT|) . we obtain the following correction to the quark 
distribution function (fTSj) : 



Af(x) 



71 



x(l-x)< 3(1 -2x) 



Li 3 (l--] -Cs 



'1 - 2x) In ( - - 1 

x 



Li 2 1- 



x 



~6 
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2 Li, 1 



x 



+ In x In ( 1 

x 
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(22) 



Notice, that 

/ Af(x) dx = (23) 
Jo 

So this correction does not change the structure function normalization 1)140. 



4 Conclusion 

We have calculated a leading-twist a s correction to the pion structure functions due to 
the final state gluon exchange between the struck quark and spectator in the Feynman gauge. 
We certainly do not claim, that it is not cancelled by the initial state interaction or does not 
vanish in some specific gauge. Although the arguments, based on the operator expansion 
of the 2-point current product are of no doubt, they are less obvious being applied to the 
4-point current correlator, considered in the sum rule approach. The spectator interaction, 
which occurs immediately before or after the virtual photon absorption, may give rise to a 
leading twist effect beyond the series (J2J). 

The initial state diagram (J^Jd) involves more scalar graphs, which require additional 
cumbersome calculations to perform. If it does not cancel the final state interaction in the 
sum, this will support the statement of jSj, and probabilistic interpretation of the structure 
functions should be reconsidered. But even if it does, the calculation of the diargams (jHt) 
with only the struck quark involved would also give a new information: its infinite part must 
reproduce the known InQ 2 contribution f)15jl. determined by the evolution, while the finite 
part would be a new correction. These diagrams altogether will provide us with refined pion 
structure functions to be fitted with experimental data. 
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Appendix 



Below we describe the calculation of the two-loop integrals (|T8j) and (|T5j). (pO)l. 

Notations and conventions. Vecotors p, g are taken to be space-like, while p + g is 
the time-like one: 

= _p 2 < o , g 2 = -Q 2 < , s = (p + g) 2 > . (Al) 

It is convenient to perform the calculation in the center of mass frame p + q = (s/s, 0). One 
may choose 

p = (vV-P 2 ,0,0,cu) , g = (v/cu 2 -Q 2 ,0,0,-cu) (A2) 

where 

/ v 2 - P 2 Q 2 

u = y , ^ = 0<?) • (A3) 

In the calculation of the 2-particle cut we use the 2-particle phase volume: 

dT (A4) 



IQnujyfs J T 
where k 2 = k\ = and 

T = -(p-k) 2 , T± = v±uj^s (A5) 



In (jA4j) we integrated out azimuthal angle <p in plane, since in the frame (jA2|) the ampli- 
tudes do not depend of it. 

The 3-particle phase volume is: 

dTs = J (2^W Wfkkl (2^f W A ^P + 1-k^k 2 -h) 
d 3 ki d 3 k 



where k\ = k\ = 0, but k 2 = A 2 , the mass A is necessary to regularize the IR divergencies. 
In spherical coordinates 



r^fs—ki pit pit p2Tr p2n 



dT 3 = j—Ys ' p ' (// ' J ' "' 7>2 ' f/ ^ J ' ,lfi - 



^2 



x5 

where we denoted 



cos a — cos 6*i cos #2 

COS(01 - 02 

sin &i sin #2 

s - 2 v / i(/c 1 + fc 2 ) + 2fcifc 2 - A 2 

cos a = — 

2kik 2 



(A7) 
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For a given 9\ the condition cos 2 
range: 

I cos a I < 1 , 



^2) < 1 sets the following restriction on the integration 



cos (6*i + a) < cos 6*2 < cos (61 — a) 
It is convenient to introduce new variable £ instead of 62 and instead of &i i2 : 

yT*, 



cos 62 = cos 61 cos a + sin 9\ sin a cos £ 



A'.; 



ti 



(A8) 



The 5-function is excluded by integrating out 02- Since the amplitudes are independent of 
the angle 0i, we integrate it also. Then the 3-particle phase volume (jA7|) takes the most 
convenient form: 



1 s r 

dT<i = -, — —r ■ — / dt] 

(2vr)4 l 6 y Aa/- 

2-box diagram consists of two cuts 

Q 



l-t L +\ 2 /s 



X 2 /(stx) 



dt 2 / sin 9i d9i I d^ 



(A9) 



V 



1 1 






A 2 


/ 

/ 


/ 




/ 




/ 1 



Cut 3 Cut 2 ( A1 °) 
and symmetric ones. Since all internal lines are massless, both cuts are IR divergent. A 
common way to regularize it is the introduction of small "photon" mass A 2 . To evaluate the 
2-particle cut, one needs to know the 1-loop box integral: 

1 



(All) 



ft I 2 (I - p) 2 (I + q) 2 [{1-P + k) 2 - A 2 ] 
The Feynman representation of the propagator product yelds the following 3-fold integral: 



16n 2 



dx / dy 



dz x 



x [Q 2 (x -y)(l-x)+ Tz(l -x) + P 2 {y - z)(l - x) - s(x - y)(y - z) + \ 2 z] (A12) 

This integral is real for s < 0, since the denominator is never zero in this case. So we 
calculate it for negative s and then analytically continuate to positive region. The result is: 



Li 2 1 



T 
P 2 



Lio 1 



T 
~Q 2 



16tt 2 sT 

Then one calculates the contribution of both 2-particle cuts: 



1 2 ST 2 7T 2 

+ 2 11 A 2 P 2 Q 2 + ~3 



(A13) 



2CuU 



dTr 



(p — k) 2 
1 

~ 128vr 3 sP 2 g 2 

4z/ , T_ 
H ^ln 



{^□| s+i0 + J nL-io} 
A 2 . , A 2 



7T" 



ln z 4 In — + 8 

s s 3 



31n^ 



PQ 



ln 5 



+ 



UJy/S 

CO 



PQ 



2 - In 



A : 



P 2 



in — m 



^ Q PQ 



p2 

Li2 [l-jT 



Y2 



Li 2 [l-|-)+^bi 



,T_ 

PQ 



1, 2 P 

+ 2 ln g 



(A14) 
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The 3-particle cut is given by the integral: 

dT 3 



Cut* 



(p + q — h) 2 (p + q — k 2 ) 2 (p - h) 2 (q - k 2 ) 2 

dT 3 (s - 2^k 1 )-\s - 2 v /jfc 2 )- 1 

[-P 2 - 2k 1 (Voo 2 -P 2 - oj cos 6 X )\ [-Q 2 - 2k 2 {sjuj 2 -Q 2 + 00 cos 6 2 )} 

Integration over the angles yields: 
1 



(A15) 



CuU 



6Att 3 s AuJy^S 

1 



1 f 1 dtt f'-^^dh 

~U 



X 



A 2 /s *1 J\ 2 /{st 1 ) 

l-t x - t 2 )P 2 Q 2 + T+ishh + P% + QH 2 ) 



(A16) 



In 



sht 2 + P 2 h + QH 2 (1 - ti - t 2 )P 2 Q 2 + T_ (sht 2 + P 2 h + Q 2 t 2 



This integral is IR divergent at £ 1)2 — > (or k 12 — > ^). The expresion in the second 
line of (|A16|) is regular there, so we put A 2 = everywhere except the integration limits; 
account of the terms neglected would lead to power corrections in A, but not logarithmic 
ones. Evaluation of (JA16j) yields exactly —1/2 of equation (|A14J) . so 



2Cut 2 + 2 Cut 3 







which proves the equation l|18)) . 
Box-triangle diagram 




Cut 3 Cut 2 

is calculated in similar way. The triangle 1-loop integral is: 



(A17) 



1 



% P (/ +p )2[(/ + A .)2_ A 2] 



16tt 2 (T-P 2 ) 
The 2-particle cut is equal to: 

i 



Li 2 1 - 



P< 



T T 
+ ln ^ ln ^ 



(A18) 



CuU 



dT, 



{k — p)'- 



i 3 I A 



2Li 3 (l-y) + 3Li 3 (l-l/y) - 2 1nyLi 2 (l-y) 



256tt 3 u^P 2 _ 
-|ln 3 y +ln^Li 2 (l-l/y) 



y+ 



T 

y = jr 2 



(A19) 
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The 3-particle cut is given by the same integral as (|A16|) . but without ti in denominator: 
Cut 3 = 



i 3 dT a 



(p + q - k 2 ) 2 (p - h) 2 (q - k 2 )' 



25Qtt 3 lu^P 2 



In In 

. y+ x . 



Li 2 (l - 1/2/h 



- 2 Li 3 (l - y+) - - In 3 y+ - (y_ <-> y 4 

In the calculation we used the following properties of the function Li3: 

Li 3 (l -z) + Li 3 (l - l/z) + U 3 (l/z) =C 3 -^\nz + l \n 3 z - 1 ln(z - 1) ln 2 z 

o 3 I 

7T 2 1 , 

H 3 (-z) - Li 3 -l/z) = - — In z - -ln 3 z 

o o 

Summing up (jA19|) and (|A20|) . one obtains the integral (fT9*|) : 
Ji = 



(A20) 



256tt 3 uj^P 2 



3 Li ; , [ 1-^V 3Li 3 fl-^ 



P- 



+ ln^Li 2 (l ^ 



7L. 



P_ 



(A21) 



The intergal (}20|) is calculated in similar way: 



J? 



P 2 



P 5 



— — - (T_ - 2T+) Li 2 1 - (T, - 2T_) Li 2 1 

- T+ In ^- In — + T_ In ^± In — 

P z s P l s 



(A22) 



In the Bjorken limit P 2 — > and 



2x x 



T_ « xP 2 . 



X 



(A23) 



Notice cancellation of the terms ~ ln"P 2 /P 2 in (|A21|) . ()A22|) : they would violate Bjorken 
scaling. 
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